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ABSTRACT 

We apply our self-consistent accretion disk corona (ADC) model, with two different geometries, to the 
broadband X-ray spectrum of the black hole candidate Cygnus X-1. As shown in a companion paper, 
models in which the Comptonizing medium is a slab surrounding the cold accretion disk cannot have a 
temperature higher than about 140 keV for optical depths greater than 0.2, resulting in spectra that are 
much softer than the observed 10 30 keV spectrum of Cyg X-1. In addition, the slab-geometry models 
predict a substantial “soft excess” at low energies, a feature not observed for Cyg X-1, and Fe Ka fluo- 
rescence lines that are stronger than observed. Previous Comptonization models in the literature have 
invoked a slab geometry with optical depth t t £ 0.3 and coronal temperature T c ~ 150 keV, but they are 
not self-consistent. Therefore, ADC models with a slab geometry are not appropriate for explaining the 
X-ray spectrum of Cyg X-1. Models with a spherical corona and an exterior disk, however, predict much 
higher self-consistent coronal temperatures than the slab-geometry models. The higher coronal tem- 
peratures are due to the lower amount of reprocessing of coronal radiation in the accretion disk, giving 
rise to a lower Compton cooling rate. Therefore, for the sphere-plus-disk geometry, the predicted spec- 
trum can be hard enough to describe the observed X-ray continuum of Cyg X-1 while predicting Fe 
fluorescence lines having an equivalent width of ~40 eV. Our best-fit parameter values for the sphere- 
plus-disk geometry are r T * 1.5 and T c & 90 keV. 

Subject headings: accretion, accretion disks — radiation mechanisms: nonthermal — radiative transfer — 
stars: coronae — X-rays: stars 


1. INTRODUCTION 

Because of its high X-ray brightness and its apparent 
large mass, Cygnus X-1 is the most studied Galactic black 
hole candidate (BHC). Although it is almost certain that the 
production of the high-energy radiation is associated with 
the accretion of matter, the exact details are uncertain. The 
source of accreting matter is most probably a combination 
of Roche lobe overflow and accretion out of the stellar wind 
of the companion, the 09.7 lab supergiant HDE 226868 
(Gies & Bolton 1986). Recent mass determinations yield a 
mass of about 18 M 0 for the companion and 10 M Q for the 
compact object (Herrero et al. 1995, in agreement with 
Dolan 1992; Sokolov 1987; Bochkarev et al. 1986; Aab et 
al. 1984; Hutchings 1978 and references therein), making 
the latter a good candidate for being a black hole. 

Cyg X-1 is usually observed in the so-called “X-ray low, 
y-ray high” state, which is roughly described by a power 
law with a photon index a ss 1.7, modified by an exponen- 
tial cutoff with an e-folding energy of about 100 keV (see 
Oda 1977; Liang & Nolan 1984; Tanaka & Lewin 1995 for 
detailed reviews of the observations). The most popular 
model used to explain the high-energy spectrum of Cyg X-1 
(and other BHCs) while in the “low state” involves inverse 
Comptonization of soft photons by a hot plasma, usually 
referred to as a corona. While the Comptonization model 
by Sunyaev & Titarchuk (1980) has been used widely in the 
past (e.g., Grebenev et al. 1993; Dobereiner et al. 1995), 
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Titarchuk (1994) showed that the typical fit parameters 
found for Cyg X-1 are outside the model’s range of validity 
(and cannot be used for interpreting the physical conditions 
in the Comptonizing medium), and he developed a more 
generalized theory. Using a composite spectrum from 
EXOSAT, Granat , and OSSE observations, Titarchuk 
(1994) found that this new accretion disk corona (ADC) 
model, having an optical depth z c % 0.6, a coronal tem- 
perature kT c ^ 150 keV, and a slab geometry, allowed for a 
rough description of Cyg X-1. 

Although Comptonization apparently accounts for the 
general behavior of the observed X-ray spectrum of Cyg 
X-1, it appears that a reprocessing component is needed for 
a full description of the spectrum. Formal descriptions of 
the different observations in terms of a reflection com- 
ponent with an underlying Comptonization spectrum yield 
better fits to the data than do pure Comptonization models. 
The presence of a weak fluorescence line from neutral iron 
(Kitamoto et al. 1990; Done et al. 1992; Marshall et al. 
1993; Ebisawa et al. 1996) and the deviation of the 2 60 
keV band from a pure power law seen by the HEAO 1 A-2 
experiment (Inoue 1989; Done et al. 1992; Gierlinski et al. 
1997), which is usually interpreted as a Compton reflection 
hump, both indicate the presence of cold or slightly ionized 
material present in the source, in addition to the Comp- 
tonizing medium. Therefore, since the radiation processes of 
Comptonization and reprocessing of coronal radiation in 
the cold accretion disk cover a large range of photon ener- 
gies ( < 2 to ~250 keV) in the predicted X-ray spectrum, it 
is necessary to use broadband spectral observations to 
place meaningful constraints on the models. 

As discussed in a companion paper (Dove, Wilms, & 
Begelman 1997a, hereafter Paper I), one serious drawback to 
solving the radiative transfer problem for an ADC with 
analytic and “linear” Monte Carlo techniques (where each 
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photon particle is propagated through a fixed background) 
is that the physical properties of the corona must be defined 
a priori. The coronal properties and the radiation field are 
strongly coupled through the processes of Compton 
cooling, Compton upscattering, and pair production and 
annihilation. Therefore, it is not clear whether the assumed 
coronal properties are self-consistent with the resulting 
radiation field. 

We have recently developed an ADC model in which the 
radiation field, the temperature and opacity of the corona, 
and the reprocessing of coronal radiation in the accretion 
disk are calculated self-consistently (Paper I). We find that 
all previous ADC models having a slab geometry and used 
to describe the X-ray spectrum of Cyg X-l are not self- 
consistent, since the coronal temperatures and opacities are 
not within the physically allowed range of combinations of 
these parameters (see Fig. 2 of Paper I). Haardt et al. (1993) 
claimed that their “ best fit ” model, having a coronal tem- 
perature kT c = 150 keV and an optical depth t t = 0.3, 
agrees within the 3 <x confidence level of their allowed kT c -x T 
relation. However, because of the quality of the data, the 3 a 
contour covers a very large region of the kT c -x T parameter 
space, of which only a very small part is within the self- 
consistently allowed region. 

In this paper, we apply our self-consistent model to Cyg 
X-l, using a more complete composite spectrum than was 
available for previous work (Wilms et al. 1996). We study 
two geometries: (1) the standard slab geometry and (2) a 
spherical corona with an exterior accretion disk. At the time 
of this writing, simultaneous X-ray spectra covering the 
range from 2 to 250 keV were not generally available. We 
therefore had to use a nonsimultaneous composite spec- 
trum from several instruments to cover this energy range 
(each observation was taken while Cyg X-l was in its “ low ” 
state). A detailed description of the individual data com- 
prised by the composite spectrum and a preliminary 
analysis are given in Wilms et al. (1996). 4 

2. THE NUMERICAL CODE 

It appears impossible to produce an analytical model of 
an ADC that properly accounts for the self-consistent 
thermal and opacity distributions of the corona while 
solving correctly the radiative transfer problem for an 
angle-dependent, high-energy radiation field in a semi- 
relativistic, nonuniform plasma (including reprocessing of 
radiation in the accretion disk). Therefore, numerical 
methods are needed. For the computations presented here, 
we use a nonlinear Monte Carlo (NLMC) code to calculate 
self-consistently the temperature and opacity of the corona, 
as well as the radiation field within and external to the 
corona. For the two geometries considered here, Compton 
scattering, photon-photon pair production, and pair anni- 
hilations are taken into account within the corona. For the 

4 We note that, because of a bug in the treatment of reprocessing of 
coronal radiation, the conclusion made in Wilms et al. (1996) that the 
slab-geometry model is able to describe the broadband spectrum of Cyg 
X-l is erroneous. As we show below, the slab-geometry model predicts a 
spectrum that is much softer than the observed spectrum. As a result of the 
bug, the yields for Compton reflection/fluorescence and thermalization 
were reversed, causing < 20% (instead of > 80%) of the coronal radiation 
to be reprocessed into thermal radiation and leading to a correspondingly 
erroneous Compton cooling rate. As discussed in Paper I, the current 
version of our code has been successfully compared to other reprocessing 
models (e.g., George, Nandra, & Fabian 1990). 


reprocessing of radiation within the accretion disk, 
Compton scattering and photoabsorption, resulting in fluo- 
rescent line emission and thermal emission, are considered. 
We enforce charge neutrality by requiring n e — n + = n p , 
where n e , n + , and n p are the number densities of electrons, 
positrons, and protons, respectively. For this paper, the 
spatial distribution of the seed opacity (the opacity not 
including the contribution due to electron-positron pairs) 
was assumed to be uniform. 

The free parameters of our models are (1) the seed opacity 
of the corona, (2) the temperature of the accretion disk as 
a function of radius, T BB (R\ and (3) the compactness param- 
eter of the corona, l c . For a set of free parameters, the 
NLMC code is iterated until the temperature of the corona, 
T c , the total opacity, x T = x p -j- 2r + , and the radiation field 
have reached a steady state. Thermal equilibrium is 
achieved when the local heating rate is balanced by the 
Compton cooling rate. The pair opacity is determined by 
requiring pair annihilations to balance pair production. 
Once the system has reached a steady state, the spectra of 
escaping radiation are recorded (in 10 bins of different incli- 
nation angles) until satisfactory statistics are obtained. A 
detailed description of the NLMC code and its use for 
ADCs is given in Paper I. 

3. SLAB GEOMETRY 
3.1. Model 

The most common geometry used to model the spectra of 
BHCs is the slab (plane-parallel) geometry. For this 
geometry, we assume that the accretion disk corona is situ- 
ated above and below an optically thick, geometrically thin, 
cold accretion disk, that all corona and disk properties are 
constant with respect to radius, and that azimuthal sym- 
metry applies. The coronal compactness parameter is 
defined as 



where z 0 is the scale height of the corona, L c is the lumi- 
nosity of the corona, R c is the radius of the corona, L Edd is 
the Eddington luminosity, R s is the Schwarzschild radius, 
is the rate of energy dissipation per unit area into the 
corona, and h = z 0 /R c . The specification of the physical 
mechanism responsible for heating the corona is not neces- 
sary. We simply assume that the energy is dissipated uni- 
formly with respect to height. For a specific value of / c , the 
self-consistent coronal properties are degenerate with 
respect to z 0 and y ¥ c (Stern et al. 1995). 

As we discuss below, the predicted amount of thermal 
excess in the spectrum of escaping radiation strongly 
depends on the temperature of the accretion disk, T BB . A 
lower limit to the temperature can be estimated by 
assuming that the disk is heated solely by illumination from 
the corona. In this case, the temperature is approximated by 
balancing absorption with thermal emission : 

°S-B ^BB “ (1 — d)F c , (2) 

where F c is the flux of coronal radiation incident onto the 
accretion disk (assumed here, by symmetry, to be equal to 
the flux of escaping radiation), <r s _ B is the Stefan-Boltzmann 
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constant, and A is the angle-averaged albedo of the disk. 
Using equation (1), the disk temperature can be expressed 
as 



( 3 ) 


We refer the reader to Paper I for a more detailed dis- 
cussion of the slab-geometry models. 

3.2. Application to Cyg X-l 

As discussed in Paper I, we find that the maximum self- 
consistent temperature of the corona is T max ~ 140 keV for 
total opacities t t > 0.2, regardless of the compactness 
parameter l c or T BB . Once this maximum temperature, T max , 
is reached, higher heating rates result in a higher opacity 
due to pair production, causing more reprocessing of 
coronal radiation in the accretion disk. Since most of the 
reprocessed radiation is emitted as thermal radiation, the 
higher amount of reprocessing gives rise to higher Compton 
cooling rates. Therefore, once pair production becomes 
important, an increase in the heating rate results in a 
decrease in the temperature. Figure 2 of Paper I shows the 
allowed parameter space of self-consistent coronal tem- 
peratures and opacities. 

In Figure 1, we show our by-eye “best fit” of the predict- 
ed angle-averaged spectrum, modified by interstellar 
absorption, to the composite spectrum of Cyg X-L For the 
slab geometry, the modeled spectrum shown is the hardest 
spectrum possible while having a cutoff energy E c > 100 


keV. It is apparent, however, that the predicted spectrum is 
still much softer than the observed spectrum. In our grid of 
models, the range of seed opacities is 0.1 < x p < 2.0, and the 
range of coronal compactness parameters is 0.1 < l c < 10 3 
(see Paper I). Out of the entire grid, the hardest spectrum 
predicted by our self-consistent models has a photon index 
of a ^ 1.8, while Cyg X-l has a power-law index of ~ 1.6 
1.7. As discussed in Wilms et al. (1996), the BBXRT 
archived data were deconvolved by assuming a power-law 
continuum, with a photon index of a = 1.62, absorbed by 
the interstellar medium with N n = 6.0 x 10 21 cm 2 . The 
COMIS/TTM data were reduced by Borkus et al. (1995), 
the HEXE data by Dobereiner et al. (1995), and the OSSE 
data by Kurfess (1995). Since the models provided reason- 
able reduced x 2 values (< 1.5), the unfolded data should be 
at least a good approximation to the actual spectrum. 
Therefore, even though we are comparing our predicted 
spectra with unfolded data, which are dependent on the 
assumed model, the result that the predicted spectrum is 
always softer than the observed spectrum is not sensitive to 
the deconvolution procedure. 

According to our linear (and therefore not self-consistent) 
Monte Carlo simulations, simulated spectra with a < 1.7, 
while T c > 100 keV, are achieved only for models with t > 
0.3. These values are consistent with past work, in which the 
spectrum of Cyg X-l has been described by ADC models 
with a slab geometry (Haardt et al. 1993; Titarchuk 1994). 
Such models, however, are not self-consistent, for the tem- 
perature and opacity values lie within the forbidden region 
of Figure 2 of Paper I. We remind the reader that the pre- 
dicted maximum temperatures are indeed upper limits, as 
our models do not include additional cooling mechanisms 
such as bremsstrahlung. (This mechanism was found to be 
negligible for the models of interest, but its inclusion would 
reduce the maximum temperatures allowed.) 

In addition to the simulated spectra having too soft a 
power-law component, our models predict a very large Fe 


10 


1 


10 ^ 


> in -1 

Cl) 10 


> 

CD 


10 L 


LU o 

^ 10' J 


10 "' 


10 


-5 ! 
0.1 



1.0 


10.0 100.0 1000.0 
E [keV] 


■i 


p IG i — slab geometry: comparison of the predicted spectrum of escaping radiation, modified by interstellar absorption, with the composite X-ray 
spectrum of Cyg X-l. For our model, T BB - 200 eV, t t = 0.28, T c = 110 keV, and N H = 6 x 10 21 cm~ 2 . 
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Ka equivalent width (EW). Our “ best fit ” slab model pre- 
dicts an EW ~ 120 eV, considerably larger than the 
observed value of EW < 30 eV (Ebisawa et al. 1996). 
Finally, the slab models also predict an excess of radiation 
for energies E < 1 keV (i.e., a “ soft excess ”). Since the “ best 
fit ” slab models are optically thin, most of the thermal emis- 
sion emitted by the disk escapes the system without inter- 
acting with the corona. We have accounted for Galactic 
absorption by using the observed column density of hydro- 
gen, N h = 6 x 10 21 cm” 2 (Wilms et al. 1996). Clearly, this 
is not sufficient to “ hide ” the soft excess for the 200 e V disk 
models. Models with a higher value of N n predict too low a 
flux for energies below 1 keV, so this soft-excess problem 
cannot be solved by having arbitrarily large hydrogen 
column densities in the model. We note, however, that the 
BBXRT observation is unique because a soft excess was not 
observed (Marshall et al. 1993). Since Cyg X-l usually has a 
soft-excess component (see, e.g., Baiucinska & Hasinger 
1991; Balucinska-Church et al. 1995; Ebisawa et al. 1996), 
the discrepancy between the predicted flux and the observed 
flux at low energies should not be considered a serious 
problem by itself. However, it does appear that our “ best 
fit ” slab-geometry model predicts a higher soft-excess com- 
ponent than observed by the ASCA GIS detector for 1993 
November 1 1 (Ebisawa et al. 1996). 

We conclude that self-consistent accretion disk corona 
models with slab geometry are not capable of reproducing the 
observed broad X-ray spectrum of Cyg X-l. 

4. SPHERE-PLUS-DISK GEOMETRY 
4.1. Model 

The main limitation of the slab-geometry models, in 
regard to explaining the high-energy spectra of Cyg X-l, is 
that there is too much reprocessing of coronal radiation in 
the cold accretion disk. Therefore, models with the repro- 
cessing matter having a smaller covering fraction (as seen 
from the corona) are more likely to explain the obser- 
vations. One such geometry is a combination of a spherical 
corona with an optically thick, geometrically thin, cold acc- 
retion disk. Here the accretion disk is assumed to be 
exterior to the corona, as shown in Figure 2. This geometry 
is very similar to the two-temperature accretion disk model 
of Shapiro & Lightman (1976), although, in our model, the 
proton temperature is assumed to be equal to the electron 
temperature. This geometry is also similar to the advection- 
dominated disk models by Narayan and collaborators 
(Narayan & Yi 1995a, 1995b; Narayan 1996, and references 
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Fig. 2. — Sphere-plus-disk geometry: spherical corona with exterior 
accretion disk. The corona has a radius R c , and the accretion disk 
(assumed to be optically thick) has an inner radius R c and an outer radius 
R d . The corona has been divided into several zones so that nonuniform 
temperature and opacity structures can be studied. For this paper, we 
assume uniform heating with volume within the corona. 


therein). However, in our models the seed photons for 
Comptonization are produced by the accretion disk 
through thermal emission rather than bremsstrahlung and 
synchrotron radiation within the corona. The accretion 
disk extends from R c to R d , where R c is the radius of the 
sphere and R d is the outer radius of the disk. We define the 
ratio of the outer disk radius to the corona radius as a = 
R d /R c . For this paper, we set a = 10, but the main results 
are not sensitive to this value (see § 4.2). 


4.1.1. Energetics 

We define the coronal compactness parameter as 

/ = — — = 4n( " A — — 
c m e c 3 R c WJIZ R c 



Similar to the slab-geometry models, we allow for intrinsic 
emission of thermal radiation by the accretion disk (in addi- 
tion to emission due to the reprocessing of coronal 
radiation). We define fP G to be the rate at which gravita- 
tional energy is dissipated directly into the corona, where 
P G is the total rate of dissipation of gravitational energy 
into the system. Consequently, (1 — f)P G is the rate of 
energy dissipation into the disk. The total luminosity of the 
accretion disk is given by 


= 0 -/)P G + L abs , (5) 

where L abs is the rate at which energy is absorbed (not 
reflected) by the disk as a result of the reprocessing of radi- 
ation emitted by the corona. Equation (5) can be expressed 
in terms of compactness parameters as 


h = ( 1 -/)/ G + / abs, ( 6 ) 

where l d = (a T /m e c 3 )LJR c , l G = ( <r T /m e c 3 )P G /R c , and / ab5 = 
( a T /m e c 3 )L ab JR c . As with the slab-geometry models, we 
constrain the models by setting (1 — f)l G = 1 for all the 
simulations. With this choice,/is given by 


As discussed in Paper I, setting (1 — f)l G to unity allows us 
to consider models where 0.01 </< 1.0, but models with 
other values of (1 — f)l G yield the same ranges of self- 
consistent coronal temperatures and opacities. 


4.1.2. Disk Temperature 

As in the slab-geometry model, we can estimate a lower 
limit to the disk temperature by assuming that the disk is 
heated solely by illumination from the corona: 

nR 2 c (a 2 - l)a s _ B <7 BB > 4 =f d (a)(l - A)LJ 2 , (8) 

where L c is the total luminosity leaving the corona and f d (a) 
is the fraction of photons leaving the corona that are repro- 
cessed in the disk. This fraction can be approximated by 
assuming that the photons leaving the corona are uniformly 
distributed over the surface of the corona with an isotropic 
distribution. With these assumptions, the fraction of escap- 
ing photons that hits the disk is given by 

v <AO> 1 f*/ 2 

m = = 2" sin (9) Amde * (9) 

J&m 
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where <AQ> is the average solid angle of the disk as seen 
from the surface of the corona. 


Afi(0) 


-f 

Jcsc 0 


2 cos 0 


a 2 -x 2 


2 + 1 - 2x sin 6 \a 2 + 1 - 2x sin 8 


dx , 


8 m - arcsin (1/a), and 6 is the angle between the normal of 
the coronal surface and the plane of the disk. In Figure 3, we 
show how the covering fraction depends on a. It is apparent 
that f d (a) is nearly constant for a > 10, and therefore the 
amount of reprocessing of radiation within the accretion 
disk is insensitive to a. The maximum value of f d is « 1/3. In 
contrast, for a slab geometry, all the downward-directed 
radiation at the base of the corona interacts with the accre- 
tion disk. Therefore the slab-geometry models predict much 
more prominent reprocessing features in the spectrum of 
escaping radiation. 

For large accretion disks (i.e., a 1), the local disk tem- 
perature can vary by more than an order of magnitude 
between the inner and the outer radius. Rather than using 
an average disk temperature for calculating the spectrum 
emitted by the disk, a more proper treatment is to take into 
account the temperature structure. We can estimate the disk 
temperature as a function of radius by equating the flux of 
absorbed coronal radiation with thermal emission, 


o's-B TUr)2nrdr = i(l - A)L c df d , (10) 


where r = R/R c , df/r)/dr is the differential covering fraction 
of a ring of radius r, and f d (r) is determined numerically by 
using equation (9). Using equation (4), the temperature can 
finally be expressed as 


t m l c \»(mR s y 

T - (r) = ,50 hrJ vbJ VrT) 


YioMq 

M 


1/4 


eV . 


(11) 


In Figure 4, we show how T BB (r) varies with r. From a 
numerical fit, we find that T BB (r) roughly decreases as T BB (r) 
oc r ~ 1 as compared with a power law of 3/4 correspond- 
ing to the standard a-disk models (Shakura & Sunyaev 
1973). It should be remembered, however, that if any accre- 
tion energy were dissipated directly into the accretion disk 
(/ < 1), then the disk temperature would be higher than the 
values given above. To address this possibility, we consider 
a simple model in which, for R > R c , the accretion disk 
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Fig. 3. — Covering fraction, of the accretion disk with an outer 
radius a (normalized to the coronal radius), averaged over the surface of 
the corona (spherical geometry). 
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temperature is calculated by assuming that the disk is heated by illumi- 
nation of coronal radiation and cooled by thermal emission (eq. [11]). 
Here l c = 25, M = 10 M 0 , and R c = 100K S . 


behaves as the standard a-disk (Shakura & Sunyaev 1973) 
while, for R < R c , all the accretion energy is dissipated into 
the corona. In order to determine the relative importance of 
the internal dissipation of energy compared with the contri- 
bution from coronal illumination, we estimate the disk tem- 
perature by neglecting the contribution from illumination. 
In addition, we make the simplifying assumption that the 
disk is optically thick and radiates as a blackbody such that 
the emitted flux is equal to the dissipation rate, 


°S-B ^BbW “ D(R) — 


/3 GMM\ 
\ 8ttR 3 ) 


I* 


( 12 ) 


where we have assumed R > R s (Frank, King, & Raine 
1992). Defining L acc = r/Mc 2 , where r\ is the efficiency by 
which the accretion process converts gravitational energy 
into radiation, and using the standard definition of the 
Eddington luminosity, L Fdd , we can express the disk tem- 
perature as 



By comparing equation (13) with equation (11) (or Fig. 4), it 
appears that the heating rate due to the direct dissipation of 
accretion energy is comparable to that from coronal illumi- 
nation. A proper treatment of the disk temperature struc- 
ture would be to equate thermal emission with both heating 
mechanisms. Such a treatment is outside the scope of this 
paper; with a heating rate equal to twice the value used in 
equation (13), the disk temperature will be roughly 1.2 times 
higher than the value given in equation (13) (T BB is pro- 
portional to the total heating rate per unit area). Therefore, 
the inclusion of direct energy dissipation into the disk 
changes the disk’s total luminosity but does not signifi- 
cantly change the temperature profile of the disk. 

For simplicity, we have assumed that the intrinsic flux of 
thermal radiation is constant with respect to radius (the 
emission due to reprocessing, however, is determined 
locally by equating emission with absorption). This approx- 
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imation is not too important, since the emission due to 
reprocessing is higher than the emission due to internal 
dissipation for the models considered in this paper. In 
regard to the soft-excess issue, for a % 10 these models over- 
estimate the predicted amount of soft excess by less than 
10 %. 


4.1.3. Covering Fraction of Corona 

In all of our models, the seed photons are produced 
within the cold accretion disk (we neglect bremsstrahlung 
radiation). For this geometry, however, not all seed photons 
have to propagate through the corona prior to escaping the 
system. Therefore, even for optically thick models, a very 
large soft excess of radiation is always predicted. The cover- 
ing fraction of the corona, evaluated on the accretion disk 
at a distance d from the center, is given by 


m) = 


AQ 

2 % 


1 r _ - r 

2 L R/R c J ’ 


(14) 


where R c is the radius of the corona. As shown in Figure 5, 
the covering fraction of the corona rapidly decreases with 
increasing R. In order to estimate the fraction of the thermal 
radiation emitted by the accretion disk that interacts with 
the corona, we average/,, over the accretion disk, weighted 
by the thermal flux, F(r): 


</c>(*) 


ji 2nrf c (r)F(r)dr 
Ji 2nrF(r)dr 


(15) 


For an accretion disk heated solely by illumination, we 
found that T(r ) oc r ~ 11 (Fig. 4). On the other hand, for the 
case in which the disk is heated solely by internal dissi- 
pation, i.e., the standard oe-disk, T(r)ocr 3/4 (Shakura & 
Sunyaev 1973). For a rough estimate of </ c )when both disk 
heating mechanisms are important, we assume 
F(r) = kT 4 (r) oc r 3 * 5 . In Figure 6, we show how this aver- 
aged covering fraction depends on a. For a > 2, we have 
<X> ^ 0.2, and at least 80% of the thermal radiation 
escapes the system without interacting with the corona, 
even if the corona is optically thick. In contrast, for the 
slab-geometry models, all the thermal radiation must pro- 
pagate through the corona prior to escaping the system. 
Therefore the sphere-plus-disk models always predict that a 
large amount of thermal radiation, relative to the amount of 
Comptonized radiation, will escape the system. Unless this 
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thermal radiation is absorbed by the interstellar matter 
(ISM), it will appear as a soft excess of radiation. 

For the sphere-plus-disk geometry, the accretion disk 
receives less illumination from the corona since the covering 
fraction is lower than in the slab-geometry case, while the 
surface area of the disk is larger. In § 4.1, we estimated the 
temperature structure due to either coronal illumination or 
internal dissipation of accretion energy (eq. [13]). For either 
case, 


T BB (r)ccr y , (16) 

where r = R/R s with y « 1.1 for heating by coronal illumi- 
nation and y = 3/4 for the standard a-disk. The flux of radi- 
ation emitted by an accretion disk is approximated by 

F{E)ccj a B^T)rdr, (17) 

where B E (T) is the Planckian distribution. In Figure 7, we 
show the corresponding spectrum, F(E\ for y — 1.1, and we 
compare F(E) to a Planckian distribution corresponding to 
the temperatures <T BB > = 60 and 80 eV. It is seen that the 
Planckian distributions approximate the high-energy tail 
fairly well. The disagreement at energies lower than 100 eV 
is irrelevant since this radiation is efficiently absorbed by 
the ISM. In addition, since the temperature decreases 
rapidly with radius, most of the radiation having energy 
E > 100 eV is emitted within the inner region of the disk, 
and the exact value of y is unimportant. Therefore, since the 
Comptonized spectrum is insensitive to the exact shape of 
the spectrum of seed photons, we used a Planckian distribu- 
tion corresponding to a single disk temperature of 80 eV to 
approximate the exact integrated disk spectrum. 

For t x > 1, the self-consistent temperature structure 
varies by roughly 30% for the case in which the heating rate 
is uniform with volume. For these optical depths, the outer 
shells of the corona are the hottest while the center region is 
the coldest. This variation is not so large that the predicted 
spectra contain a “hardening 1 ’ feature that mimics the 
hardening due to Compton reflection (Skibo & Dermer 
1995). 


4.2. Application to Cyg X-l 

We computed a grid of models in which the range of seed 
opacities is 0.5 < x p < 4.0 and the range of coronal com- 
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Fig. 7. — Spectrum of thermal radiation emitted from an accretion disk 
with a radial temperature structure, T BB (r), as given by eq. (11) ( solid line). 
Here l c = 23, R c = 100K S , and M - 7 M Q . The dashed line is a Planckian 
distribution corresponding to a temperature T BB = 60 eV; the dash-dotted 
line is a Planckian distribution corresponding to a temperature T BB = 80 
eV. All distributions have been normalized to unity. 


pactness parameters is 0.1 < l c < 100. In Figure 8, for 
/cT bb = 50 eV, we compare the predicted spectrum of escap- 
ing radiation, modified by interstellar absorption, of our 
“best fit ” model to the composite spectrum of Cyg X-t. 
Since the seed photons enter the corona from the exterior of 
the sphere, the effective optical depth of the sphere is 
reduced. Therefore the total optical depth [defined by t T = 
R c <r T (n e + n p \ where n e = n p + n+] must be higher than 
that for the slab geometries in order for the models to 
predict similar power laws. Our “ best fit ” model has a total 
optical depth r x = 1.5 and an average coronal temperature 
<T C > = 90 keV. 


With this geometry, the corona is able to reach much 
higher temperatures as compared with the slab models, 
because the corona is “photon starved/’ i.e., the luminosity 
of the seed photons is much less than the luminosity of the 
hard X-rays (Zdziarski, Coppi, & Lamb 1990). As discussed 
above, most photons that are reprocessed within the accre- 
tion disk do not reenter the corona. Therefore the Compton 
cooling mechanism that prohibits slab-geometry coronae 
from having high temperatures is not as efficient in keeping 
the spherical coronae cool. In Figure 9, we show the 
allowed regime of optical depths and average coronal tem- 
peratures for the sphere-plus-disk geometry. Since the 
coronae are allowed to be much hotter than the slab- 
geometry models, there are self-consistent (7^., re- 
combinations such that the observed power law and cutoff 
in BHCs can be reproduced. Since more than 80% of the 
thermal radiation escapes the system without passing 
through the corona, the soft excess is very large for these 
models. As in the case with the slab-geometry models, 
Galactic absorption is not sufficient to “hide” the soft 
excess for the 200 eV disk models. In Figure 8, the predicted 
thermal excess appears to be consistent with the observ- 
ational data, and the entire broadband spectrum is ade- 
quately described by the model. As discussed in the 
Appendix, even though we are fitting “ unfolded ” data, we 
are confident that this model provides a good representa- 
tion of the intrinsic photon spectrum of Cyg X-l. However, 
the main point of this paper is to motivate a more detailed 
study of the sphere-plus-disk model, using a more rigorous 
method of data analysis. 

A nice feature of the sphere-plus-disk model is that it 
naturally predicts an “effective” disk temperature to be 
much lower than in the slab-geometry models. In fact, 
because of the lower disk temperature, this model appears 
to predict a soft excess that is consistent with the obser- 
vations (Balucinska & Hasinger 1991; Balucinska-Church 
et al. 1995). The quality of the observations, however, is not 
too good for E < 2 keV, which is unfortunate since this is 
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Fig. 8— Comparison of the predicted spectrum of escaping radiation, modified by interstellar absorption, to the X-ray spectrum of Cyg X-l. This is for 
the sphere-plus-disk geometry. Here T BB = 50 eV, t t = 1.5, T c — 90 keV, and N H — 6 x 10 21 cm 2 . 
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Fig. 9. — Allowed temperature and opacity regime for self-consistent 
ADC models with a sphere-plus-disk geometry. The solid line is derived 
from a “ fit by eye ” to the numerical results. For a given total optical depth, 
temperatures above the solid line are not possible. For all models, the 
blackbody temperature of the disk is kT BB = 50 eV. Different symbols 
represent models with different seed opacities, while different points having 
the same symbol represent different coronal compactness parameters. 


the energy range in which the spectrum is most sensitive to 
the disk temperature. Also, as noted in § 3.2, Cyg X-l 
usually contains a soft-excess component. We plan to 
compare this model with better observational data in this 
energy range (e.g., Rossi X-Ray Timing Explorer [ RXTE ] 
and ASCA data) in order to constrain this model further. 

To test whether the predicted thermal radiation from the 
accretion disk is observable in the ultraviolet, we have com- 
pared the ultraviolet flux predicted by our model for the 
spectral band from 1200 to 1950 A with the flux observed in 
this band by IVE. Our predicted flux is 2-3 orders of mag- 
nitude smaller than the flux in an WE low-dispersion spec- 
trum of HDE 226868 made in 1980 June. This is consistent 
with earlier findings that the ultraviolet spectrum of the 
HDE 226868-Cyg X-l system is dominated by emission 
from the O star (Treves et al. 1980; Pravdo et al. 1980). 
Therefore, WE observations cannot be used to constrain 
the disk temperature of our models. 

Because of the smaller covering fraction of the accretion 
disk, the models predict weaker reprocessing features in the 
escaping radiation field as compared with the slab models. 
In principle, in modeling BHCs, the value of a can be con- 
strained by the equivalent width of the Fe Ka fluorescence 
line, which is proportional to the solid angle of the disk, 
f d (a). For Cyg X-l, the EW of the Fe Ka line is measured to 
be EW < 70 eV (Ebisawa et al. 1996; Gierlinski et al. 1997). 
With a = 10 (corresponding to a disk covering fraction f d % 
0.3), our best-fit models predict an EW of ~60 eV. There- 
fore models with smaller f d will predict EWs that are too 
low to be consistent with the observations. For f d > 0.25, we 
have a > 3 (Fig. 3). On the other hand, since fj^a) does not 
vary significantly for a 10, this method cannot provide an 
upper limit on the size of the accretion disk. 

5. DISCUSSION 

Our self-consistent slab accretion disk corona models are 
unable to explain the broadband X-ray spectra of Cyg X-l. 
The modeled coronae are either too cold or have an optical 
depth that is too small, resulting in spectra that are much 


softer than the observed spectrum. All previous ADC 
models that have successfully described the spectra of Cyg 
X-l have temperature and opacity values that are outside 
the allowed region of self-consistent values. In addition, as 
discussed in Paper I, the predicted angle-averaged equiva- 
lent width of the Fe Ka fluorescence line is EW > 150 eV, a 
value roughly 3 times higher than the measured value for 
Cyg X-l (Ebisawa et al. 1996), indicating that the covering 
fraction of the reprocessing material is Eijln ~ 0.3. 

We believe that there is no way around these shortcomings 
of the slab-geometry model and that these ADC models are 
not appropriate for explaining the high-energy spectra of 
BHCs. This claim is in agreement with Gierlinski et al. 
(1997), who, by fitting the joint Gm^a-OSSE observation of 
Cyg X-l with a power law plus a reflection component, find 
that the solid angle of the cold medium is ~0.3(2rc), signifi- 
cantly lower than the angle corresponding to the slab 
geometry, a result consistent with Ebisawa et al. (1996). 
Gierlinski et al. (1997) also argue that the corona must be 
“ soft-photon starved,” which rules out a slab geometry. 

For ADC models with a sphere-plus-disk geometry, there 
are self-consistent temperature and opacity combinations 
such that the power-law and cutoff portions of the spectra 
of BHCs can be explained. In addition, as a result of the 
small covering fraction of the disk as seen from the corona, 
it is possible that the “effective” disk temperature (a flux- 
weighted average disk temperature) can be as low as < 80 
eV. With such low values, the models predict a thermal 
excess (or lack thereof) that appears to be consistent with 
the observations of Cyg X-l, since Galactic absorption is 
more efficient at these lower energies. It is still unclear, 
however, whether an ADC model having this geometry is 
stable, as the Shapiro & Lightman (1976) model suffers from 
a thermal instability (Piran 1978). In our models, we can 
only assume a corona/disk morphology, but we cannot 
determine whether this morphology is thermally or dynami- 
cally stable. The sphere-plus-disk geometry does approx- 
imate the geometry corresponding to the advection- 
dominated models of Narayan (1996, and references 
therein), which appear to be thermally stable. In addition, 
our best-fit models for Cyg X-l have an optical depth t t % 
1.5 and a coronal temperature kT c % 90 keV, values that are 
consistent with the advection-dominated models if the mass 
accretion rate M ~ (0.1-0.3)M Edd , where M Edd is the 
Eddington accretion rate (Narayan 1996). 

The results obtained here for Cyg X-l are also applicable 
to other low-state BHCs. In Figure 10, we compare the 
spectra of the black hole X-ray transients GS 2023 + 338 
and GRO J0422 + 32 to the predicted spectra from our 
sphere-plus-disk geometry model. Occasionally, black hole 
X-ray transients have spectra that are dominated by a 
strong soft excess at lower energies and a pure power-law 
component at high energies. In the low state, however, GS 
2023 + 338 ( = V404 Cyg) and GRO J0422 + 32 did not 
contain a soft excess, and their spectra are very similar to 
that of the low state of Cyg X-l. Therefore it appears that 
photon-starved ADC models, such as the sphere-plus-disk 
geometry, are applicable to explaining the observed X-ray 
spectra of many BHCs. 

Our “fits” to the unfolded data should be interpreted as 
a motivation to consider this model in much more detail. 
For example, we have recently implemented our sphere- 
plus-disk model in tabular form for the XSPEC user- 
defined model. We have recently implemented the 
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Fig. 10. — Comparisons between the composite spectra of GRO J0422 4- 32 (Kroeger et al. 1996) and V404 Cyg (in ’t Zand et al. 1992; Dobereiner et al. 
1994; 1989 June 20 spectrum) and the predicted spectra from the best-fit sphere-plus-disk models. The spectra have been scaled arbitrarily. For V404 Cyg, the 
parameter values of our best-fit model are kT BB = 52 keV, r T = 3.4, and N H = 2.4 x 10 22 cm 2 . For GRO J0422 + 32, kT BB = 75 keV, r T = 1.9, and 
N h = 1.6 x 10 22 cm 2 . 


sphere-plus-disk geometry spectral grid as a user-defined 
model for XSPEC (version 9; Arnaud 1996), and we are 
currently analyzing new data on Cyg X-l and other BHCs 
from the RXTE satellite. Our preliminary results are in 
agreement with the results presented in this paper (Dove et 
al. 1997b). 
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APPENDIX 

Comparisons of our models with observations have relied on the photon spectrum presented in Figure 8. Although the use of 
“ unfolded ” X-ray spectra is very common in the astronomical literature, it is important to keep in mind how these spectra are 
prepared. Because of the poor energy resolution of today’s X-ray detectors, the inversion of the detector response matrix is 
generally impossible (Blisset & Cruise 1979 and references therein). Thus, in order to obtain an estimate of the intrinsic 
photon spectrum from a source, a spectral model has to be fitted to the spectrum observed by the X-ray detector, usually 
using a x 2 reduction method and an estimate for the detector response matrix (Gorenstein, Gursky, & Garmire 1968). If the x 2 
value that results from the fitting process is small, then it is assumed that the modeled photon spectrum is a good representa- 
tion of the observed photon spectrum. In order to produce an “unfolded” spectrum, the additional assumption is made that 
the residual (for a given detector energy channel) between the modeled spectrum (folded through the response matrix) and the 
observed count rate is a good representation of the deviation between the spectral model and the intrinsic photon spectrum 
(in energy space). Thus, at a given energy, the product of the modeled photon flux and the ratio between the observed count 
rate and the (estimated) model count rate is what is called the “ unfolded ” spectrum. 

Although a small x 2 value of the fit to the “ unfolded ” spectrum usually implies that the model is an adequate description of 
the intrinsic spectrum, sometimes this is not true. Fitting our models to the observed data in detector space (using the method 
described above) would, in principle, be the best method of determining the “goodness of fit.” However, the CPU time needed 
to obtain model spectra with a signal-to-noise ratio that is suitable for spectral fitting is prohibitively high, and we were forced 
to use a simplified approach; for a particular detector, we simulated an observation using our modeled spectrum as the 
intrinsic spectrum, the correct detector response matrix, background estimate, and observational Poisson noise. This simu- 
lated spectrum was then fitted with several spectral models that are typically used to describe the X-ray spectra of BHCs (e.g., 
a power law with an exponential cutoff). 

We performed these tests using several typical X-ray detectors ( ROSAT , RXTE PCA, HEXE) and using typical count rates 
and observation times for these detectors. In all the tests, the fit parameters found by using our model spectrum were in 
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agreement with typical fit parameters found in the real observation of Cyg X-l. In the simulated RXTE PCA observation, 
e.g., the model spectrum below 30 keV can be roughly described by an absorbed power law with an index of a « 1.62. This 
index is similar to those found, e.g., by Done et al. (1992) and Dobereiner et al. (1995). In the residuals of the simulation, the Fe 
line and the reflection hump are clearly visible. Using a more appropriate model that incorporates reflection, the simulated 
spectrum is best described by reflection of a power law of index T % 1.72 off a cold accretion disk with covering factor 0.4. 
Using an assumed length of 5 ks for the ROSAT PSPC observation, it is possible to see evidence for the soft excess in the 
simulated observation. According to this simulated observation, the soft excess as seen by the ROSAT PSPC can be described 
by a blackbody with a temperature of about 100 eV or less, which is roughly comparable to the observations by Balucinska- 
Church et al. (1995). Therefore we are confident that the model spectrum of Figure 8 is a good representation of the intrinsic 
photon spectrum of Cyg X-l. 
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